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Abstract —  Small-scale  variation  in  temperature  and  salinity 
can  lead  to  localized  changes  in  the  index  of  refraction  and  can 
distort  electro-optical  (EO)  signal  transmission  in  ocean  and 
atmosphere.  This  phenomenon  is  well-studied  in  the  atmosphere 
and  in  this  context  is  generally  called  “optical  turbulence”.  Less 
is  known  about  how  turbulent  fluctuations  in  the  ocean  distort 
EO  signal  transmissions,  an  effect  that  can  impact  various 
underwater  applications,  from  diver  visibility  to  active  and 
passive  remote  sensing.  To  provide  a  test  bed  for  the  study  of  the 
impacts  from  turbulent  flows  on  EO  signal  transmission,  as  well 
as  to  examine  and  mitigate  turbulence  effects,  we  set  up  a 
laboratory  turbulence  environment  allowing  the  variation  of 
turbulence  intensity.  Convective  turbulence  is  generated  in  a 
large  Rayleigh-Benard  type  tank  (5m  by  0.5m  by  0.5m)  and  the 
turbulent  flow  is  quantified  using  a  suite  of  sensors  that  includes 
high-resolution  Acoustic  Doppler  Velocimeter  profilers  (Vectrino 
Profiler)  and  fast  thermistor  probes  (PME  Conductivity- 
Temperature  probe).  These  measurements  allow  the 
characterization  of  turbulent  kinetic  energy  and  temperature 
variance  dissipation  rates  in  the  tank,  for  different  convective 
strengths.  Optical  image  degradation  in  the  tank  is  then  assessed 
in  relation  to  turbulence  intensity.  The  turbulence  measurements 
are  further  complemented  by  very  high-resolution  computational 
fluid  dynamics  simulations  of  convective  turbulence  emulating 
the  tank  environment.  These  numerical  simulations  supplement 
the  sparse  laboratory  measurements,  providing  full  fields  of 
temperature  and  velocity  in  the  tank.  The  numerical  data 
compared  well  to  the  laboratory  data  and  both  conformed  to  the 
Kolmogorov  spectrum  of  turbulence  and  the  Batchelor  spectrum 
of  temperature  fluctuations.  The  numerical  model  was  able  to 
qualitatively  reproduce  the  turbulence  fields  observed  in  the 
laboratory  tank.  Quantitatively,  the  numerical  simulations  are 
consistent  with  the  observed  £  in  the  tank,  but  do  not  fully  resolve 
the  temperature  gradients  and  thus  underestimate  /.  The  unique 
approach  of  integrating  optical  techniques,  turbulence 
measurements  and  numerical  simulations  can  help  advance  our 
understanding  of  how  to  mitigate  the  effects  of  turbulence 
impacts  on  underwater  optical  signal  transmission,  as  well  as  on 
the  use  of  optical  techniques  to  probe  oceanic  processes. 

Keywords — optical  turbulence;  turbulence  measurements; 
ADV;  Rayleigh-Benard  tank;  turbulent  dissipation  rate 


I.  Introduction 

The  term  “optical  turbulence”  is  widely  used  in 
atmospheric  optics  to  describe  the  phenomenon  of  distortion  of 
EO  transmission  due  to  variations  in  the  index  of  refraction 
caused  by  changes  of  air  temperature  along  the  optical  path. 
This  is  similar  to  the  visual  perturbations  that  can  be  observed 
over  a  hot  road  or  burning  candle.  A  similar  phenomenon 
occurs  underwater,  and  here,  can  be  caused  by  variations  in 
either  temperature  or  salinity  that  are  associated  with  turbulent 
microstructure.  Most  often,  temperature  fluctuations  are  the 
dominating  factor  affecting  the  index  of  refraction  [1],  except 
in  the  case  of  strong  freshwater  or  salt  water  influence,  such  as 
in  river  outflows  or  estuaries,  or  possibly  in  surface  lenses 
generated  by  rainfall.  In  this  study,  we  neglect  the  influence  of 
salinity  and  instead  focus  on  the  effect  of  temperature. 

Underwater  “optical  turbulence”  at  sea  was  investigated  in 
the  1970s  by  [2],  who  used  this  term  to  describe  “small 
inhomogeneities  in  the  index  of  refraction  of  seawater,  their 
origins,  and  the  effects  they  have  on  underwater  optical 
systems”.  The  phenomenon  has  since  received  limited 
attention,  despite  its  potential  to  affect  a  wide  range  of 
applications,  from  diver  visibility  to  active  and  passive  remote 
sensing  of  the  ocean.  More  recently,  [1]  looked  at  the  effect  of 
“light  scattering  on  oceanic  turbulence”  with  numerical  studies 
compared  to  measurements  in  a  small  laboratory  tank.  Two 
recent  field  studies  aimed  at  characterizing  naturally-occurring 
“optical  turbulence”  in  the  aquatic  environment  highlight  the 
difficulties  associated  with  collecting  concurrent  data  on  optics 
and  turbulence  in  the  ocean  [3]  or  lakes  [4]. 

The  impact  of  turbulent  fluctuations  on  optical  signal 
transmission  can  be  quantified  using  the  so-called  optical 
turbulence  parameter,  Sn,  which  is  the  oceanic  equivalent  of  the 
atmospheric  optical  turbulence  coefficient  Cn2  [5]  [6].  Sn  is  a 
function  of  turbulent  kinetic  energy  and  temperature  variance 
dissipation  rates,  8  and  %,  respectively: 

sn  ~  x e  1/3  (1) 
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Thus,  in  order  to  estimate  Sn  and  thus  the  amount  of  image 
degradation  from  turbulence,  we  need  information  on  both 
velocity  and  temperature  fluctuations.  These  measurements 
need  to  be  of  sufficient  resolution  to  allow  inferring  dissipation 
rates  8  and  %. 

To  simplify  the  measurements  and  reduce  the  number  of 
parameters  involved,  such  as,  f.  ex.,  platform  motion  polluting 
the  velocity  data  and  scattering  due  to  particles  degrading  the 
optical  signal,  which  affect  field  data  in  particular,  we 
developed  a  laboratory  environment  for  the  study  of 
underwater  “optical  turbulence”.  The  laboratory  setup  also 
allows  us  to  vary  turbulence  intensity  and  thus  generate  a 
“controlled  turbulence  environment”  in  terms  of  parameter 
space.  The  laboratory  experiments  are  supplemented  by 
numerical  simulations  using  computational  fluid  dynamics 
(CFD).  The  numerical  simulations  provide  full  fields  of 
temperature  and  velocity  and  thus  provide  a  better  view  of  the 
large-scale  flow  field  and  distribution  of  turbulence  parameters 
than  could  be  gathered  with  the  sparse  laboratory 
measurements  alone.  In  this  paper,  we  present  the  setup  and 
results  from  laboratory  and  numerical  experiments  in  this 
turbulence  environment  designed  for  the  study  of  underwater 
“optical  turbulence”. 

II.  METHODS 

A.  Laboratory  Setup 

The  laboratory  setup  consists  of  a  large  acrylic  tank  which 
is  5m  long  and  has  a  cross  section  of  0.5m  by  0.5m  and  is 
outfitted  with  stainless  steel  plates  at  the  bottom  and  top  that 
can  be  temperature  controlled  (Fig.  1,  left).  In  this  tank, 
convective  Rayleigh-Benard  type  turbulence  is  generated  by 
heating  and  cooling  the  bottom  and  top,  respectively.  The 
strength  of  the  convective  turbulence  in  the  tank  is  a  function 
of  the  temperature  difference  across  the  tank  and  can  be 
characterized  in  terms  of  the  Rayleigh  number,  defined  as 

Ra  =  gaATd3/(vDT)  (2) 

Where  g  is  the  acceleration  due  to  gravity,  a  is  the  thermal 
expansion  coefficient,  AT  is  the  temperature  difference 
between  the  plates,  d  is  the  distance  between  the  plates,  v  is  the 
kinematic  viscosity,  and  DT  is  the  thermal  diffusivity. 


Fig.  1.  Laboratory  tank  at  NRLSSC.  The  photo  on  the  right  shows  the 
Vectrino  profiler  ADV  and  the  CT  temperature  probe. 
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Fig.  2.  Example  of  the  “Numerical  Tank"  used  to  simulate  Rayleigh-Benard 
convection  and  emulating  the  laboratory  tank  setup.  The  temperature  field  (in 
K)  is  shown  and  several  convective  plumes  are  visible. 

The  Rayleigh  number  of  the  flow  can  be  changed  by 
changing  the  plate  temperatures,  and  thus  the  turbulence 
intensity  can  be  varied.  In  our  experiments,  Ra  ranges  from 
1.5- 10 10  to  around  4T010,  corresponding  to  a  temperature 
difference  AT  between  the  plates  of  6K  and  16K,  respectively. 

The  turbulence  in  the  tank  is  quantified  by  high-resolution 
Acoustic  Doppler  Velocimeter  profilers  (Nortek  Vectrino 
Profiler)  and  fast  thermistor  probes  (PME  high-resolution 
conductivity-temperature  (CT)  probe)  (Fig.  1,  right).  These 
instruments  provide  high-resolution  velocity  and  temperature 
measurements,  at  100  and  64Hz,  respectively.  Three  AD  Vs 
and  two  CT  probes  were  mounted  in  the  tank  and  collected 
time  series  of  high-resolution  velocity  and 
temperature/conductivity  for  the  subsequent  estimation  of  s 
and  %.  Data  were  collected  at  a  sampling  frequency  of  100Hz 
with  the  Vectrino  Profiler  and  at  64Hz  with  the  CT  probes, 
which  were  controlled  by  a  Nortek  Vector  ADV  in  our  setup. 
Turbulent  kinetic  energy  dissipation  rate  s  and  temperature 
variance  dissipation  rate  x  were  calculated  from  the  velocity 
and  temperature  measurements  via  spectral  fitting  to 
Kolmogorov  spectra  (for  velocity)  and  Batchelor  spectra  (for 
temperature)  and  compared  to  values  obtained  from  the 
numerical  simulations  of  convective  turbulence  in  the  tank  for 
comparable  Rayleigh  number  and  setup  [7]  [8]  [9]  [10].  The 
turbulence  data  was  put  into  the  context  of  measurements  of 
optical  target  clarity,  by  placing  a  high-speed  imaging  camera 
and  active  optical  target,  an  iPad  displaying  optical  resolution 
charts,  at  opposite  ends  of  the  tank,  providing  an  optical  path 
length  of  around  5m.  To  quantify  the  extent  of  image 
degradation  from  optical  turbulence,  an  image  quality  metric, 
namely  the  Structural  Similarity  Index  Method  (SSIM)  was 
applied  [11].  The  SSIM  measures  the  similarity  between  two 
images,  where  one  is  considered  to  be  of  “perfect  quality”. 

B.  Numerical  Tank 

The  sparse  laboratory  measurements  (three  Vectrino 
Profilers  and  two  temperature  probes)  were  complemented  by 
CFD  simulations  of  the  convective  tank.  These  three- 
dimensional,  very  high-resolution,  non-hydro  static  numerical 
simulations  provide  full  fields  of  temperature  and  velocity  for 
the  estimation  of  turbulence  parameters  and  their  impact  on  the 
optics. 


This  project  was  supported  by  ONR/NRL  73-4951/73-6604.  Silvia  Matt 
was  supported  by  a  NRC  Research  Associateship  and  a  NRL  Karle 
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Fig.  3.  Model  results  at  time  t  =  1050s,  showing  a  sideview  of  temperature  T 
(top),  vertical  velocity  (Uz)  (middle)  and  the  Index  of  Refraction  (IOR) 
calculated  from  the  temperature  field  (bottom).  The  data  is  shown  on  the 
center  plane  of  the  tank. 
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Fig.  4.  Model  results  at  time  t  =  1050s,  showing  streamlines  in  the  domain 
colored  by  vertical  velocity  Uz  (top,  in  m/s).  The  bottom  plot  shows  the  same 
but  looking  at  it  from  the  other  side  of  the  tank  (x-  and  y-axis  reversed).  The 
updrafts  and  downdrafts  illustrate  that  there  is  a  circulation  in  the  tank  in  the 
cross-sectional  direction. 
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Fig.  5.  Image  degradation  for  different  “optical  turbulence”  strengths.  The 
images  are  taken  with  an  optical  path  of  5m  (the  length  of  the  laboratory  tank) 
of  an  optical  resolution  chart  across  a  region  of  turbulence.  Left  is  for  a  still 
tank  (no  turbulence),  the  middle  is  for  a  “strong”  level  of  turbulence,  and  on 
the  right  is  the  highest  level  of  turbulence  (“extreme”)  we  can  achieve  in  the 
tank. 


The  numerical  experiments  were  performed  with  the  open- 
source  CFD  package  OpenFOAM  using  a  Large-Eddy 
Simulation  (LES)  approach.  In  LES,  the  larger-scale  eddies  in 
the  flow  are  explicitly  resolved,  while  the  scales  smaller  than 
the  grid-size  are  modeled  [12].  The  traditional  Smagorinsky 
model  was  chosen  as  the  sub-grid  scale  model  [13].  Here,  we 
present  results  from  a  very  high-resolution,  millimeter-scale 
simulation  with  Ax  =  Ay  =  5mm,  Az  =  2.5mm,  which 
corresponds  to  20  million  grid  points  in  the  5m  by  0.5m  by 
0.5m  domain  (Fig.  2).  Exploratory  simulations  at  lower 
resolution  (Ax  =  Ay  =  Az  =  1  cm)  were  run  on  a  modern,  high- 
performance,  dual  six-core  Linux  desktop,  whereas  the 
production  runs  at  the  millimeter-scale  resolution  required 
High-Performance  Computing  resources  at  the  DoD 
Supercomputing  Resource  Center,  due  to  the  high 
computational  cost. 

III.  Results 

Temperature  and  velocity  from  the  numerical  simulations 
provide  a  view  of  the  overall  circulation  in  the  tank,  thus 
supplementing  the  sparse  laboratory  measurements.  The 
simulations  illustrate  the  convective  cells  that  are  established 
in  the  tank  (Fig.  3).  The  size  and  number  of  the  convective 
cells  that  develop  in  the  tank  are  a  function  of  the  tank 
dimensions,  in  particular  the  tank  height,  since  the  water  rises 
and  sinks  and  gets  diverted  once  it  reaches  the  solid 
boundaries  at  the  top  and  bottom.  With  a  tank  of  depth  d  = 
0.5m  and  length  L  =  5m,  we  observe  on  the  order  of  ten 
convective  cells  in  our  domain.  In  addition  to  the  convective 
cells,  secondary  circulations,  namely  in  the  cross-sectional 
direction,  develop  in  the  domain  (Fig.  4).  These  circulations 
can  be  visually  confirmed  in  the  laboratory  when  adding  a 
tracer  to  the  fully  developed  flow  field,  such  as,  for  example, 
the  seeding  material  needed  to  collect  ADV  data. 

The  model  temperature  fields  and  the  index  of  refraction 
(IOR)  calculated  from  these  fields,  using  the  empirical  relation 
described  in  [14],  provide  an  illustration  of  the  disturbance 
expected  to  be  encountered  by  an  optical  beam  passing 
through  the  tank  (Fig.  3). 


Images  from  the  laboratory  show  the  image  degradation 
due  to  optical  turbulence  for  a  case  of  “strong”  (Fig.  5,  center, 
AT  ~  6K)  and  “extreme”  optical  turbulence  (Fig.  5,  right,  AT 
~  16K).  Here,  particle  scattering  is  secondary  to  the  changes 
in  the  index  of  refraction  due  to  temperature  micro  structure. 
Note  that  the  effect  of  optical  turbulence  is  more  pronounced 
at  the  higher  spatial  frequencies.  When  applying  the  SSIM 
metric  to  a  sequence  of  video  images  in  order  to  quantify  the 
extent  of  image  degradation  from  optical  turbulence,  we 
confirm  the  differing  amounts  of  image  degradation  for  the 
two  turbulence  cases.  The  time  series  also  hints  at  a  difference 
in  the  frequency  of  variations  between  the  two  cases. 


SSIM  image  degradation  metric 


Fig.  6.  Image  degradation  for  different  (optical)  turbulence  strength  as  shown 
in  Fig.  5  quantified  by  applying  the  SSIM  metric. 

To  quantify  the  turbulence  in  the  tank  and  to  estimate  the 
turbulence  parameters  necessary  to  estimate  the  optical 
turbulence  coefficient  Sn,  we  calculate  turbulent  kinetic  energy 
dissipation  rates  8  and  temperature  variance  dissipation  rates  % 
from  the  data  collected  in  the  laboratory.  These  values  are 
then  compared  to  8  and  %  calculated  from  the  numerical 
experiments.  The  numerical  data  compared  well  to  the 
laboratory  data  and  both  conformed  to  the  Kolmogorov 
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Fig.  7.  Energy  spectra  (left)  and  temperature  gradient  spectra  (right)  from 
laboratory  and  model  data.  The  energy  spectra  are  shown  only  for  the  highest 
level  of  turbulence  (compare  to  Fig.  5,  right),  since  the  spectra  are  not 
significantly  changed  for  different  convective  turbulence  strengths.  The 
temperature  gradient  spectra  do  resolve  the  difference  in  turbulence  strengths. 

spectrum  of  turbulence  and  the  Batchelor  spectrum  of 
temperature  gradients  (Fig.  7).  The  numerical  model  was  able 
to  qualitatively  reproduce  the  turbulence  fields  observed  in  the 
laboratory  tank.  Quantitatively,  the  numerical  simulations  are 
consistent  with  the  observed  8  in  the  tank,  despite  the  fact  that 
they  do  not  resolve  the  spectrum  down  to  the  Kolmogorov 
microscale,  which  is  on  the  order  of  mm  for  this  flow. 
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Fig.  8.  Turbulent  kinetic  energy  dissipation  rate  8  from  velocity  section  along 
the  length  of  the  tank  at  each  x-z-point,  in  W/kg.  Left  row:  low  Ra  turbulence 
at  times  t  =  250s  (top)  and  t  =  1675s  (bottom).  Right  row:  high  Ra  turbulence 
at  times  t  =  250s  (top)  and  t  =  1675s  (bottom). 


Extreme  turbulence  -  Sensor  P 


The  laboratory  data  show  that  even  for  convective  strength 
with  dramatically  different  impact  on  the  optics  (Figs.  5,  6),  the 
s  in  the  tank  remains  within  one  order  of  magnitude.  To 
illustrate  how  s  is  expected  to  vary  across  the  tank  cross- 
section,  we  can  take  the  spectrum  from  the  numerical  model 
data  at  every  5m-long  velocity  section  in  the  tank,  for  each  x- 
and  z-position  (Fig.  8). 
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The  value  of  s  stays  mostly  within  one  order  of  magnitude 
across  the  tank.  The  variation  near  the  boundaries  due  to 
boundary  layer  effects  can  be  expected,  and  is  consistent  with 
the  cross-sectional  circulation  seen  in  the  velocity  field.  While 
the  s  from  the  model  is  consistent  with  the  laboratory  value, 
both  are  of  O(10"7  W/kg)  for  the  experiments  shown,  the 
numerical  simulations  do  not  fully  resolve  the  temperature 
gradients  and  thus  underestimate  %,  even  at  the  high  resolution 
used  in  our  experiments  (A  x  =  A  y  =  5mm,  A  z  =  2.5mm;  20 
M  grid  points).  The  values  of  %  from  both  the  laboratory  and 
the  numerical  simulations  show  a  significant  difference  for 
different  convective  turbulence  strengths,  but  since  the 
estimates  from  the  numerical  simulations  are  strongly 
resolution  dependent  and  lower  than  the  laboratory  values  for 
comparable  turbulence  strength,  estimating  the  model  sub-grid 
scale  contribution  to  %  is  of  critical  importance  and  is  the 
subject  of  ongoing  work.  Since  the  effect  on  the  optics  is 
driven  by  changes  in  the  index  of  refraction  due  to  temperature 
variations,  these  results  reemphasize  the  importance  of 
characterizing  in  detail  the  temperature  distribution  to  assess 
the  impact  of  the  turbulent  fluctuations  on  the  optics. 

IV.  Conclusions 

In  order  to  study  the  impact  of  temperature 
microstructure  on  underwater  optical  signal  transmission,  we 
performed  experiments  in  a  controlled  laboratory  environment 
complemented  by  high-resolution,  non-hydrostatic  numerical 
simulations.  The  goal  was  to  develop  a  setup  where  turbulence 
levels  can  be  controlled  and  fully  characterized.  This  setup, 
which  allows  for  repeatable  experiments  under  controlled 
conditions,  can  help  us  understand  processes  involved  in 
optical  turbulence  and  provide  a  platform  for  the  testing  of 
optical  techniques  to  mitigate  turbulence  effects  underwater. 
Optical  turbulence  is  mainly  due  to  temperature  (or  salinity) 
variations  affecting  the  IOR  of  water,  and  to  adequately 
describe  the  effect  on  the  optics,  it  is  particularly  critical  to 
resolve  the  temperature  gradients.  This  can  present  a  challenge 
in  both  the  laboratory  and  the  model,  due  to  noise  and 
resolution  requirements,  respectively.  Further  work  is  needed 
to  address  questions  related  to  sub-grid  scale  contributions  in 
LES  to  the  rate  of  temperature  variance  dissipation.  The 
unique  approach  of  integrating  optical  techniques,  turbulence 
measurements  and  numerical  simulations  can  help  advance 
our  understanding  of  how  to  mitigate  the  effects  of  turbulence 
impacts  on  underwater  optical  signal  transmission,  as  well  as 
on  the  use  of  optical  techniques  to  probe  oceanic  processes. 
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